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Abstract 

The 3-d inverse scattering problem of the reconstruction of the unknown dielec¬ 
tric permittivity in the generalized Helmholtz equation is considered. The main 
difference with the conventional inverse scattering problems is that only the mod¬ 
ulus of the scattering wave field is measured. The phase is not measured. The 
initializing wave field is the incident plane wave. On the other hand, in the previous 
recent works of the authors about the “phaseless topic” the case of the point source 
was considered [201EI1E2]. Two reconstruction procedures are developed for a lin¬ 
earized case. However, the linearization is not the Born approximation. This means 
that, unlike the Born approximation, our linearization does not break down when 
the frequency tends to the infinity. Applications are in imaging of nanostructures 
and biological cells. 
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1 Introduction 

This paper is the continuation of three recent publications [201 ED E2] of the authors 
dedicated to the reconstruction procedures for 3-d Phaseless Inverse Scattering Problems 
(PISPs). Our goal here is to extend the result of [22] to the case when the wave propaga¬ 
tion process is generated by an incident plane wave rather than by the point source of izg. 
We end up our derivation with the inversion of the 2-d Radon transform. This can be done 
by the well known technique, see, e.g. [2B]. We also propose a different inversion method 
based on the integral geometry. The idea of the second method goes back to the famous 
paper of Allan Cormack published in 1963 [6]. It is well known that his publication [6] 
led him to the Nobel prize in 1979, along with Godfrey Hounsheld, for his work on X-ray 
computed tomography (CT), see https://en.wikipedia.org/wiki/Allan_McLeod_Cormack. 
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In 120] a reconstruction procedure was proposed for a 3-d PISP for the Schrodinger 
equation in the frequency domain. The next question is about the reconstruction for the 
generalized Helmholtz equation. An important difference between these two is that in the 
generalized Helmholtz equation the unknown coefficient is multiplied where A; > 0 is 
the frequency, and k is varied in our problems. On the other hand, in the Schrodinger 
equation the unknown potential is not multiplied by k"^. This makes the reconstruction 
procedure for the PISP for the Schrodinger equation easier than that for the generalized 
Helmholtz equation. 

In [21] the reconstruction procedure was developed for the case when the Born approx¬ 
imation is considered for the generalized Helmholtz equation. However linearization of 
the Born approximation breaks down for large values of k. Hence, in [22] the reconstruc¬ 
tion procedure for that equation was developed for the case when the wave propagation 
process is generated by the point source. Even though a linearization was used in iza, 
this linearization does not break down for large k. Still, we point out that both here and 
in iza the phase of the scattered wave field is reconstructed without the linearization, 
separately for each direction of the incident wave field. The only approximation we use 
when reconstructing the phase is that we ignore the term O {1/k) for k —>■ cxo. 

Inverse scattering problems without the phase information arise in imaging of struc¬ 
tures whose sizes are of the micron range or less. Recall that 1 micron (1/rm) = 

For example, nano structures typically have sizes of hundreds on nanometers, which is 
about 0.1/rm. Therefore, the wavelength for this imaging should be about 0.1/im. The 
second example is in imaging of living biological cells. Sizes of cells are between b/im and 
100/im [531131]- Either optical of X-ray radiation is used in this imaging. It is well known 
that for the micron range of wavelengths only the intensity of the scattered wave field can 
be measured, and the phase cannot be measured [ZllHl El EH HQ]. The intensity is the 
square modulus of the scattered complex valued wave field. We point out that standard 
statements of inverse scattering problems assume that both the modulus and the phase 
of the complex valued scattered wave field can be measured, see, e.g. [mi El Ea ESI Eg. 

As to the uniqueness results PISPs, see mEiEi] for the 1-d case and HZIEIE] 
for the 3-d case. However, proofs of these results are not constructive. The first rigorous 
reconstruction procedures were proposed recently by the authors [221 El ES] and indepen¬ 
dently and for different statements of problems by Novikov [2D1IST] . In [T21E] numerical 
methods for reconstructions of shapes of obstacles from phaseless scattered data were 
developed. As to some heuristic reconstruction algorithms, we refer to, e.g. works of 
physicists [ZllEl El 1321 HQ] • 

In section 2 we state both forward and inverse problems. In section 3 we study 
geodesics lines generated by the refractive index n{x),x G In section 4 we study 

an auxiliary Cauchy problem for a hyperbolic equation. In section 5 we show how to 
reconstruct the phase for high frequencies. In section 6 we linearize the problem. Finally 
in section 7 we present two reconstruction procedures. 

2 Statements of forward and inverse problems 

Let R > 0 and H = {|x| < B} C be the ball of the radius B with the center at {x = 0}. 
Denote S = {|a;| = R} . Let the refractive index n{x),x G be a function satisfying the 
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following conditions 


n{x) G n{x) = 1 + (5{x), (2.1) 

I3{x) > 0, I3{x) = 0 for X G \ 11. (2.2) 

The smoothness of n{x) will be explained below (see the proof of Theorem 1). Let 
= {z/ G : |z/| = 1} be the unit sphere. We consider the following equation 

Am + fc^?7,^(x)M = 0, X G (2.3) 

u (x, k) is the complex valued wave field and fc > 0 is a positive real number. We seek the 
solution of equation fl2.3|) in the form 

u (x, k, p) = exp {—ikx ■ p) + Usc {x, fc, p) , (2.4) 

where the function Ugc {x, k, p) satisfies the radiation condition at the inhnity, 

du 

—^ + ikUsc = o{r~^), r = |x| —)■ cx). (2.5) 

or 

It is well known that the problem fl2.3l) - fl2.5p has unique solution u (x, k, p) G (M^), Vz/ G 
5^ [5]. Furthermore, theorem 6.17 of [H] implies that, given fl2.ip . the function u (x, k, p) G 
^ 16 +a g X), Vz/ G 5^, where (M^) is the Holder space. 

We model the propagation of the electric wave field in by the solution of the problem 
fl2.3p - fl2.5p . This modeling was numerically justified in |1] in the case of the time domain. 
It was shown numerically in |1] that this modeling can replace the modeling via the full 
Maxwell’s system, provided that only a single component of the electric field is incident 
upon the medium. Then this component dominates two others and its propagation is well 
governed by the single equation fl2.3p . This conclusion was verihed via accurate imaging 
using electromagnetic experimental data in, e.g. Chapter 5 of [3] and unsg. 

We consider the following inverse problem: 

Phaseless Inverse Scattering Problem. Consider a number k^ > 0. Find the 
function (3 (x) assuming that the following function f (x, p, k) is given 

f (x, z/, k) = \usc {x, z/, fc) 1^ , Vx G S', Vz/ G S*^, Vfc > k^, (2.6) 


3 Geodesic lines 

Our reconstruction procedure is essentially using geodesic lines generated by the function 
?7.(x).Hence, we consider these lines in this section. The Riemannian metric generated by 
the function n(x) is 

dr = n{x) \dx \, \dx\ = \/ (dxi)^ + {dx2Y + {dx^y. (3.1) 

For each vector p E 3“^ define the plane F(z/) as 

E(z/) = {e G : e ■ z/ = -B}. (3.2) 

Let .^Q G F(z/) be an arbitrary point. Then — ^q) ■ p = 0, V.^ G Fl(z/). Hence, the plane 
E(z/) is orthogonal to the vector p. Also, note that 

S(z/)nO = 0. (3.3) 
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We can represent an arbitrary point .^g G S(i/) as 


^0 = V 2 ) = -By + + 62772 , (^ 1 , V 2 ) e 

where nnit vectors u, Ci, 62 form an orthogonal triple. 

Let the fnnction v) be the solntion of the Canchy problem for the eikonal eqnation, 


|V5r(^,z/)| = n(0, r|sH = 0, 


snch that 




> 0 if • 17 > —5, 
< 0 if ■ 17 < -B. 


(3.4) 


(3.6) 


It is well known that |r(.^, v) \ is the Riemannian distance between the point ^ and the 
plane Physically, |r(.^, Z 7 )| is the travel time between ^ and the plane Ll(z 7 ). For 

^ ■ u < —B the fnnction v) has the form = ^ ■ v + B. To hnd this fnnction for 

^ • Z 7 > —B, we need to solve the problem 03.41) in this domain. It is well known that to 
solve this problem, one needs to solve a system of ordinary differential eqnations. These 
eqnations also define geodesic lines of the Riemannian metric. They are (see, e.g. EH): 


- P “^P _ V7 (I lew _ 1 

ds nmy da da 


(3.6) 


where s is a parameter and p = V^r(^, u). Consider an arbitrary point ^ 0 ( 771 , 772 ) ^ 
and the solntion of the eqnations 03.61) with the Canchy data 


= ^o(’,l.P2). Pl—0 = n(^o(Pl.P.2))P. t|,=O = 0. 


(3.7) 


Here the parameter s > 0. The solntion of the problem 03.61) . 03.7p defines the geodesic 
line which passes trongh the point ^ 0 ( 771 , 772 ) in the direction u. Hence, this line intersects 
the plane Fl(z 7 ) orthogonally. For s > 0 the solntion determines the geodesic line ^ = 
g{s, 77 ;^, 772 , y) and the vector p = h{s, 772 , y) which lays in the tangent direction to the 
geodesic line. It is well known from the theory of Ordinary Differential Eqnations that if 
the fnnction n{x) G C'™'(M^), m > 2, then fnnctions g, h are C"* ^—smooth fnnctions of 
their variables. 

By (1321) 




(977 1 


s=0 


= ei. 




dV2 


s=0 


= 62 - 


Moreover, it follows from fl3.6p and fl3.7p that 


ds 


s=0 


Hence, 


ni^oiVi,V2)) 

9(^1, ^2: ^3 


= V. 


^(5,771,772) 


= 1 ^ 0 . 


(3.8) 


s=0 


By fl3.8p the eqnality x = g{s, 77 ^^, 772 , y) can be nniqnely solved with respect to s, 77 ^^, 772 
for those points x which are snfficiently close to the plane S(i/), as 

s = s{x,v), 77 ^ = 77 i((r,z/), 772 = 772 ( 37 , 7 /). 
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Hence, the equation 


^ = 9{s, Vi{x, V 2 i^^ 9{s, x,u), s e [0, s{x, u)] 

defines the geodesic line r(a;, u) that passes through points x and ^Q^rj^^x, i'),rj 2 {x, u)) = 
^q{x, v). Extend r(a;, v) ioi ^ ■ u < —B by the equation ^ = .^0(^5 + ^^1 s < 0. This 

line intersects the plane S(z/) at the point ^ orthogonally. The Riemannian 

distance between points x and ^)- Note that functions g{s,x, u) 

and h{s, x, u) = h{s, ri^{x, u),r] 2 {x, z/), u) are smooth functions of their arguments. 

Since, moreover, h{s{x,h'),x,h') = Vxt{x,i'), we conclude that t(x, z/) is C™'-smooth 
function. In our case t{x, v) is C^^-smooth function of x for x enough closed to S(z/). 

We have constructed geodesic lines r(a;, z/) above only “locally”, i.e. only for those 
points X which are located sufficiently close to the plane E(z/). However, we need to 
consider these lines “globally”. Hence, everywhere below we rely on the following As¬ 
sumption: 

Assumption. We assume that geodesic lines of the metric h3.1\) with conditions ^3.4\ ), 
Lt. ,51) satisfy the regularity condition in M^. In other words, for each vector z/ G and 
for each point a: G there exists a single geodesic line T (x, z/) connecting x with the 
plane E(z/) such that T (a:, z/) intersects S(z/) orthogonally. 

It is well known from the Hadamard-Cartan theorem [2] that in any simply connected 
complete manifold with a non positive curvature each two points can be connected by a 
single geodesic line. The manifold (hi, n) is called the manifold of a non positive curvature, 
if the section curvatures K{x, ex) < 0 for all a; G fl and for all two-dimensional planes a. 
A sufficient condition for K{x, ex) < 0 was derived in [39] 

t > 0 . V. e a ^ e (3.9) 

In our case if condition fl3.9p holds for x G fl, then is also holds for all x G M^. 

Our Assumption implies that for any z/ G 5^ and for any x G there exists a unique 
pair (,^q(x, z/), s(x, z/)) with ,^q(x, z/) G E(z/) such that the geodesic line r(x, z/) connects 
points X and intersects the plane E(z/) orthogonally. Hence, there exists a 

function g{s, x, v) such that the equation of the geodesic line r(x, z/) is given by 

r(x, z/) = : i = g{s, X, z/), s G [0, s(x, z/)]} 

and X = g{s{x, z^), x, z/). Then s(x, z/) = r(x, u) is the Riemannian length of r(x, u). Also, 
p(s(x, z^), X, u) = Vxt{x, z/) and |Va;r(x, z/)| = rz,(x). 

4 Auxiliary Cauchy problem for a hyperbolic equa¬ 
tion 

Modifying ideas of [201 EH Eg, we consider in this section an auxiliary Cauchy problem 
for a hyperbolic equation with the incident plane wave. Consider the hyperbolic equation 

n^{x)vu — An = 0, x G f G M. (4.1) 
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We seek the solution v {x, t, u) of equation fl4.ip in the form 

V = 6{t — X ■ u) +v {x, t, u ), (4.2) 

where 6{t—x-u) is the incident plane wave propagating in the direction u, and the function 
V {x, t, v) satisfies the following initial conditions 

V \t<-B= 0. (4.3) 


Let T > 0 be an arbitrary number. Let (p{x, v) be the solution of the following problem 


V^ip{x,vf\ = n^{x), 


(4.4) 


(p(a:, v) = X ■ u for all x ■ v < —B. 

Let the function A{x, v) be 


(4.5) 


A(x, v) 


exp I n 2(^)A^(p(^, u)dT , x-iy> -B, 

r{x,u) ) 

1, X ■ V < —B. 


(4.6) 


The reason of the second line of fl4.6p is fl4.5p as well as the fact that r(x, u) is the straight 
line for x ■ v < —B. Let T > 0 be an arbitrary number. Consider the domain D{i>,T) 
defined as 

D{iy, T) = {(x, f)| max(—5, ip{x, u)) < t <T}. (4.7) 

Let H (t) be the Heaviside function, 


H{t) 


l,t> 0, 

0,t < 0. 


Theorem 1. Let conditions as well as Assumption hold. Let z/ G 

he an arbitrary vector and T be an arbitrary number. Then there exists unique solution 
v{x,t,J^) of the problem which can he represented in the form 


v{x, t, v) = A{x, v)5{t — (p(x, v)) + n(x, t, v)H{t — (p(x, z/)), x G t G (—C) 0 , T), (4.8) 


where the function H(x, z/) is defined in ( [^.6] ), the function n(x, t, z/) G C‘^{D{h',T)) and 
v{x, t, z/) = 0 for t < —B. 

Proof. Introduce functions 6k{t) as 


-/•rC 

e_,it) = 5" it), e.^it) = 5' it), e_^{t) = 5{t), e^it) = -Hit), fc = o, i, 2 ,... 


Observe that 6*).(t) = 9k_iit) for all k > —2. We seek the solution to the problem fl4.ip - 
fl4.3p in the form 


n(x, t, u) = H(x, iy)0_i (t — (p(x, u)) + E afe(x, z/)6'fc(t - (fix, z/)) + Wrix, t, z/) (4.9) 

k=0 
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for X G G M,r > 1. Substituting representation fl4.9p in fl4.ip . using \V = 
n^{x) and equating coefficients at dk{t), we obtain formulas for finding coefficients A{x, u), 
ak{x, z/): 


2VA{x, u) ■ V(p(x, u) + A{x, z/) A(p(x, u) = 0, 

2Vao{x, v) ■ V(p(x, v) + ao{x, z/)A(p(a;, v) = AA{x, z/), (4-10) 

2Vak{x, u) ■ V(p{x, u) + ak{x, h')Aip{x, u) = Aak-i{x, u), k > 1. 

Since by fl4.2p and 04.31) v {x,t,^) = S{t — x ■ u) for t < —B, then we obtain 

A{x, u) = 1, ior X ■ u < —B, (4-11) 

(p{x, u) = X ■ u, for X ■ v < —B, (4-12) 

ak{x, u) = 0, k > 0, for x ■ u < —B. (4-13) 

Moreover using 04.1l) - 04.3p and 04.9I) - 04.11|) . we obtain the following Cauchy problem 
for the residual Ws{x, t, u) of the expansion 04.9p 

n^{x)d^Wr — AWr = Fr{x, t, z/), (4-14) 

Wr\t<-B = 0, (4.15) 

Fr{x, t, v) = {Aar{x, v))9r{t — (p{x, v)). (4.16) 

We now construct solutions of equations 04.inp with the Cauchy data 04.1ip - 04.13p . 
We rewrite equations 03.61) of the geodesic line r(x, z/) with the vector p(x, z/) = V(^(a;, z/) 
in the following form 

dx _ 9 , , , , dp(x,u) „ , ,, 

— = n ^{x)p{x,iy), - - -= V(lnn(a;)), (4-17) 

dr dr 

where r is the arc Riemannian length with the element length given by the formula 
dr = n(^)|(i^|. Hence, we have along the geodesic line r(a;, u) 

2'VA{x, u) ■ V(p{x, ly) = 2VA{x, v) ■ p{x, v) = 2n'^{x)'VA{x, z/) • ^ = 2n^{x) ^^^^'' . 

dr dr 

Hence, the hrst equation 04.101) becomes 

2n^(a;) ^^^f’ + A{x, iy)A(p{x, u) = 0. (4.18) 

dr 

The solution of this equation with the initial condition 04.lip is given by 04.6p . 

Then equations for ak{x, z/), k > 0, can be rewritten in the form 


d / ak{x,iy) \ _ Aafc-i(^, z/) 
dT\A{x,iy)) 2n‘^{x)A{x,i')^ 


(4.19) 


where we should formally have a_i(x, z/) = A{x, u). It is easy to verify 04.19p . if expressing 
A(p{x, z/) via A{x, z/) using 04.18p . 
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Taking into acconnt the initial data fl4.13p . we obtain 


ak{x, u) = A{x, u) 




Aafc i(^, I/) 


dr, fc > 0, 


X ■ u > —B. 


(4.20) 


Let m > 1 be a snfficiently large integer which we will choose below. If n{x) G 
then we have 

^{x,iy) e C'™(M3),kl(x,z/) e C^-\R^),ak{x,u) e C^-^-^\R^),ar{x,u) E C^-^-^^{R^). 

Hence, the function f, z/) G C\D{i>,T)), where I = niin(m — 6 — 2r, r). By 04.13p 
and 04.16P 

Fr{x, t,iy) = 0 for t < ip{x, u) and for x • z/ < —B. (4.21) 

We now prove that 

Wr {x, t,u) = 0 for t G [—B, ip {x, z/)] . (4.22) 

Indeed, 04.15P implies that 04.22p is true for t < —B. Let now (x^,t^) be an arbitrary 
point of the domain Q(iy) = {(x,t) ^ R'^ '■ —B < t < (p(x, z/)}. Let r(x,x°) be the 
Riemannian distance between points x and x°. Denote 

/L(z/, x°) = {(x, t) G — B < t < (p(x°, v) — r(x, x°)}. 

Hence, iL(z/, x°) is the inner part of the characteristic cone with the vertex at the point 
(x°, (p(x°, z/)) and bounded by the plane t = —B. Obviously (x°, t°) G K{u, x°). By 04.211) 
Fr{x,t,^) = 0 for (x,f) G iL(z/, x°). Apply now the method of energy estimates to the 
problem 04.141) . 04.15p in the domain K{i>, x°). Observe that the resulting surface integral 
will be non-negative on the lateral boundary t = (p{x^, v) — r(x, x°) of the characteristic 
cone K{v,x^). Also, the integral over {t = —B} O K{u,x^) will be zero due to 04.15p . 
Hence, we obtain Wr{x,t,^) = 0 for all (x,f) G K{i/,x^). Thus, Wr{x^,t^,i^) = 0. Since 
{x^,F) is an arbitrary point of the domain Q(z/), then we conclude that 04.221) holds. 

Furthermore, using theorems 3.2, 4.1, corollary 4.2 and energy estimates of Chapter 4 
of [23], one can easily prove that there exists unique solution Wr of equation 04.141) . 04.15|) 
such that 

Wr{x, —B, u) = dtWr{x, —B, z/) = 0 

and Wr{x, t, v) G {t, u, T)), dtWr{-, t, v) G H^{Y{t, u, T)) for all t G {—B, T), where 

Y{t, u, T) = D{u, T) n {t = const}. Choosing m = 15 and r = 3, we obtain I = 3. Thus, 
Ws{x,t,^) e H'^{Y{t,u,T)) and dtw^{x,t-,v) G H^{Y{t,u,T)). Therefore the embedding 
theorem implies that W3 G (D(z/, T)). 

Setting 

, (t — 09(x, z/))^ 

v{x, t,n) = ^ ak{x, v) --h tz;3(x, t, u), 

k=0 

we obtain 04. 8 p as well as the required smoothness v{x,t, G C^(D(z/, T)). □ 




























4.1 Connection with the problem (l2.3l) -( l2T5l) 


Let v{x,t, u) be the solution of the problem fl4.ip - fl4.3p which is guaranteed by Theorem 
1. Fix an arbitrary bounded domain $ C And consider the behavior of the functions 
d^v{x, = 0,1,2 for a: G <F and t — )■ oo. We refer to Theorem 4 of Chapter 10 of the 

book |13] as well as to Remark 3 after that theorem. It follows from these results that 
these functions decay exponentially as f ^ oo as long as a: G <h. 

Consider the Fourier transform V {x, k, y) of the function v. 


V{x,k,v)= / V {x,t, u) exp (—ikt) dt, x G G. 


(4.23) 


Next, theorem 3.3 of |13] and theorem 6 of Chapter 9 of H guarantee that V{x, k, v) = 
u{x, k, v), where the function u{x, k, u) is the solution to the equation fl2.3^ - 112.4p . 

Using the representation fl4.8p . we obtain 

OO 

u{x,k,h')=A{x,h')exp{—ikip{x,h'))+ j v {x,t,h') exp {—ikt) dt, x G <F. (4.24) 

ip{x,u) 

Denote 

div+{x,(p{x,iy),iy) = lim div{x,t,iy),j = 0,1. 

Integrating by parts in fl4.24p . we obtain for x G $ 

/ , \ .i/ X / •, / exp(-ikip(x,iy))^ , , ^ ^ 

uix, k, v) = A(x, v) exp(—*fc(p(x, u)) H- - - i'), v) 

ik 


+ — / Vt(x,t,iy)exp(—ikt)dt. 

ik J 


Hence, 


u{x, k, u) = A{x, v) exp(—ifc(p(x, z/)) + O ^ — ) , A; —)• oo, Vx G <F, Vz/ G S'^. 

This and 02.41) imply that for all x G D, z/ G 

Usc{x, k, u) = A{x, y) exp{ikip{x, z/)) — exp(iA;x • z/) + O ^ — ) , k ^ oo. (4.25) 
For every vector u E denote {u) = {x E S : x ■ u > 0}. 


5 Approximate reconstruction of the function Usc{x, k, u) 
for all pairs z/ G 5^, x G *S'^(u) 

Ignoring the term O {1/k) in 04.25^ and using 02.6p . we obtain the following approximate 
formula 

/(x, k, u) = A^(x, u) + 1 — 2A{x, u) cos [k{^p{x, z/) — x • z/)], Vz/ G S'^, Vx G (5.1) 
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We now fix the pair v e 5^, x G and consider f{x, k, y) as the function of k for 

fc > fci, where ki > /cq is a sufficiently large number. It is possible to hgure out whether 
or not (p(x, z/) = x ■ u. Indeed, it follows from fl5.ip that (p{x, u) = x ■ u if and only if 
f{x, k, p) = const, for k > ki, i.e. if and only if dkf{x, k, p) = 0, V/c > ki. 

Suppose now that f{x, fc, p) ^ const, for k > ki. It follows from flS.ip that there exists 
a number k 2 > ki such that 

f{x, k 2 , p) = max/(x, k, p) = {A{x, p) + 1)^ . 

k>ki 

In particular, we hnd from here the number A{x, p), 

A{x,p) = y'f{x,k 2 ,p) - 1. (5.2) 

By fl5.ll) there exists a sequence {kn {x, p)}'^^^ C (^2 {x, p) , oo) such that 
f{x, kn {x, p) ,p) = maxf{x, /c, p) and k 2 < k^ < ...kn < ... 

k>ki 


Hence, 

ks (x, p) {(p{x, p) — X ■ p) = k 2 (x, p) {(p{x, p) — X ■ p) + 271. 
Thus, we approximate the number (p(x, z/), 

27r 


(p(x, p) = X ■ P + 


ks (x, p) - k 2 (x, p) 


(6.3) 


Thus, it follows from fl4.25p that formulae 05.21) and 05.31) provide us with an ap¬ 
proximation of the function Usc{x, k, p) for sufficiently large values of k and for all pairs 
p G S'^,x G S^{p). We have obtained this approximation only using the data 02.61) for 
our inverse problem. We use the word “approximation” here because we got 05.1|) via 
ignoring the term O {1/k) in 04.25p . Next, we should reconstruct the function (3 (x). This 
is done in section 7. 


6 Linearization 

We assume in sections 6 and 7 that 


ll/^llc(!rz) << (6-1) 

Use the linearization method for (p{x, p) proposed in [21| (see chapter 3 in this book and 
mm)- Then, we represent solution to problem (??) in the form 

(p{x,p) = (Pq{x,p) + ip.i^{x,p), ( 6 . 2 ) 

where (Pq{x, p) = x-p. By 02. ip n^ (x) = 1 -|- 2(5 (x) -1- (5“^ (x). By 06. 1|) the term (3“^ (x) can 
be neglected in the latter expression. Substituting the representation 06.2p in (??) and 
neglecting terms (V(pi)^ and we hnd 


V 2 :(pi(x, p) ■ Va;(Po(^) = /^(^)) p) = 0 for all x ■ p < —B. (6.3) 
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Since Vx^q{,x, u) = u, then the left hand side of equation fl6.3p is the derivative in the 
direction u. Hence, 

0 

= j f3{x + va)da, (6.4) 

-x-u—B 

The right hand side of 06.41) is the integral over the segment / (x, z/) of a straight line and 
a is its arc length. The segment I (x, z/) is stretched from the point x in the direction —u 
and continues until reaching the point x (x, z/) = x — (x • z^ + H) z/ G E (z/). Note that by 
03.3p X (x, u) ^ H. 

Let u E and x G *S'’''(z/). Denote hy y = |/(x, z/) = x — 2(x ■ z/)z/ the second 
intersection point of the straight line ^ = x + ua with S. Let L (x, z/) denotes the segment 
of that straight line connecting points x and y (x, u). Since the function (x) = 0 outside 
of the ball D and since S = dfl, then equation 06.4p is equivalent with 

(P]^(x, z/)= J (3 {^) da,Wi' E S‘^,\/x E (6.5) 

L{x,u) 


7 Reconstructions 

Formula 06.5p enables us to use two reconstruction procedures for finding the function 
/3 (x). The first one is the inverse 2-d Radon transform, and the second one solves integral 
equations of the Abel kind. As it was pointed out in Introduction, the idea of second 
procedure goes back to the paper [ 6 ]. In both cases we can reconstruct the function (3 (x) 
separately in each 2 -d cross section {xs = a = const.} of the ball D. 

For any number a G M consider the plane Pa = {xa = a}. Consider the disk Qa = 
D n Pq and let the circle Sa = S D Pa he its boundary. Then the radius of this circle 
is Ba = VB'^ — a^. Clearly 7 ^ 0 for a G {—B,B) and Qa = 0 for |a| > B. Denote 
ha = (0, 0, a) G Qa the orthogonal projection of the origin on the plane Pa- We have 

B B 

D= IJ Qa,dn:=S= IJ 

a=—B a=—B 

Denote 

Sq = = (i^i,«^2,i^3) e : z/3 = 0} , 

^ai^) = {x E Sa ■ X ■ U > 0}, Vz/ G Sq. 

7.1 Reconstruction using the inverse Radon transform 

In this subsection we present the reconstruction formula based on the inverse 2-d Radon 
transform. First, we parameterize L (x, z/) in the conventional way of the parametrization 
of the Radon transform [20] • For z/ G Sq,x E let m be the unit normal vector 

to the line L (x, z/) lying in the plane Pa and pointing outside of the point Oa. Since we 
work now only with the plane Pa, we discard the third component. Let a G (0, 27r] be the 
angle between m and the xi—axis. Then m = m (a) = (cos a, sin a). Let d be the signed 
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distance between L {x, u) and the point 0^ (page 9 of [2H])- It is clear that there exists a 
one-to-one correspondence between pairs {x, u) and (m (a) ,d),d (x, u) G {—Ba, Ba ), 

(x, z/) (m (a) ,d)]i> E S'g, x G a = a (x, z^) G (0, 27r] , d = d (x, z/). (7.1) 

Hence, we can write 


L (x, z/) = = {zi, Z 2 , a) eQa- z- m{a) = d} , (7.2) 

where = ( 2 : 1 , Z2) G and parameters a = a{x, u) and d = d{x, u) are dehned as in 

dZU). 

Consider an arbitrary fnnction q = q{z) G (Pa) snch that q {z) = 0 for 2 : G Pa\Qa- 
Hence, 

J q{z)da= j q (z) da, Wx E Sa{iz), E Sq, (7.3) 

L{x,i') z-m(a)=d 

where a = a (x, z/), s = s (x, z/) are as in fl7.ll) . In fl7.3l) a is the arc length and the 
parametrization of L{x,iy) is given in fl7.2l) . Therefore, nsing fl7.1l) - fl7.3p . we can dehne 
the 2-d Radon transform Rq of the fnnction q as 

(Rq) (x, u) = (Rq) {a, d) = j g (z) da. (7.4) 

z-m{a)=d 


Let R~^ be the operator which is the inverse to the operator R in fl7.4p . The specihc 
form of R~^ is well known, see, e.g. [22|. Hence, we do not present this form here for 
brevity. Then formnlae fl6.5l) and fl7.4p imply that 

/3 {z,a) = R~^ z/)) {z,a) , Va G {-B,B) , V 2 ; G Qa- (7.5) 

By fl5.3|) and fl6.2l) the fnnction 9 ?^(x, v) is approximately known for all v E and for 
all X G 5'+(z/). Thns, formnla fl7.5p completes onr reconstrnction process via the inverse 
Radon transform. 

7.2 Reconstruction via solutions of integral equations of the 
Abel type 

Consider now again an arbitrary nnmber a G {—B,B) and an arbitrary pair u G *S'q,x G 
^^(z/). Then the point ?/ (x, z/) = {yi,y 2 ,a) E Sa- Given the above pair x, z/, there exists 
nniqne point y (x, u) E Sa- And vice versa: given a point y E Sa and a vector 1 / E S^, there 
exists nniqne point x G *S'+(z/) snch that the segment L (x, v) passes throngh points x and 
y. Hence, we denote below M (x, y) the segment of the straight line passing throngh points 
x,y E Sa- Clearly the set {M [x,y)}^y^g^ coincides with the set {L 
i.e. {M {x,y)]^ y^g^ = {L (x, . Keeping this in mind, we rewrite ([S3]) as 

'il^{x,y)= j /3 {^)da,\/x,y E Sa,\/a E {-B,B), (7.6) 

M{x,y) 
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where the function ijj{x, y) is constructed from the function (pi{x, v) in an obvious manner, 
i.e. for each pair x,y & Sa we find the vector v & such that x G 5+ (z/) and y = y {x,^), 
and then set 'ip{x, y) = U ^)) = ^)- 

In the plane Pa we introduce polar coordinates r,(p of the variable point ^ = (^ 1 ,^ 2 ) 
as = r cos 0,^2 = rsinc/). We characterize M{x,y) by the polar coordinates (p, a) of 
its middle point {x + y)/2. Hence, \x — y\ = Change variables in the 

integral fl7.6p as 

a r,a = \/B'^ — a? — — y/r^ — p^. 

Then da = —rdr/— P^- The equation of M(x, y) can be rewritten as 

0 = a + (—1)-^ arccos-, j = l,2, r > p, (7.7) 

r 

where j = 1 corresponds to the part of the segment M(x, y) between points (x + p)/2 and 
X, and j = 2 to the rest of M(x, y). 

Obviously, there exists a one-to-one correspondence, up to the symmetry mapping 
(x,p) -v^ iy^x) between pairs x,y E Sa and (p, a) G (0, i?) x (0,27r). Denote 

/9(r COS0, r sin 0, a) = /3{r, p, a) and 'ip{p, a, a) = 'ip(^x, y). 

Using fl7.7p . we rewrite equation fl7.6p as 

9 Ba 

^ ^ /* — ■ p r — 

V / /j(r,a-I-(-1)^ arccos-,a) -.dr ='i(:(^p,a, z). (7.8) 

^7 r - p^ 

•’ p 

Represent functions /3(r, p, a) and V’(p, 0 ) via Fourier series, 

00 

^(r,p,a)= ^ ^„(r,a)exp(m0). 


ij{p,a,z)= V’n(P,a)exp(ma). 

n=—oo 

Multiplying both sides of fl7.8p by exp(— ma)/(27r) and integrating with respect to a G 
(0, 27r), we obtain for all n = 0, ±1, ±2,... 

Ba 

/ '^ / p\ vdv 

/3„(r,a)cos (^narccos-j --^=== = V’n(p, a), P e (0, Ra) (7.9) 

p 

This is the integral equation of the Abel type. To solve equation 07.91) . we apply hrst the 
operator A to both sides of 07.91) . where 


A(h (p)) (cn) 


1 f Kp)pdp 

Tri ^p2 -a;2’ 

S 


U G (0, Ba) ■ 
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Changing the limits of the integration, we obtain 


r 

— , ^, COS f narccos -]dp 

v/p2 - a;2 . ^^2 _ p2 V r) 

U) 

= ^ (V’n(P,a)) M- 

Change variables in the inner integral fl7.10p as 

p ^ 9,p^ = oj"^ cos^ {9/2) + sin^ {9/2). 



LJ 


dr 


Then 


2pdp = (r^ — sin 9 cos 9d9, 

■\/p2 — a;2 . ^^2 _ p2 ^ Q Q 

Hence, eqnation flT.lOp can be rewritten as 


Ba 


/3„(r, a)Qn{r, u)dr = 2A ( i)^{p, a)) {u ), 


(7.10) 


(7.11) 


Qn{r,uj) 



n arccos■ 


\/r2 cos2(6*/2) + ci;2sin^(6*/2) 


d9. 


We have Qn{uj,u;) = 1. Hence, differentiating fl7.1ip with respect to u, we obtain Volterra 
integral eqnation of the second kind 


Ba 


d r. 


a) - / /^n(c a)Tn{r, Uj)dr = ^2A (V'„(p, ^)) (c^) , a; G (0, H„), (7.12) 


Tn{r,uj) = 


n 


7i\/r^ — 

cos^(6'/2) + a;^sin^(6'/2) 


:X 


TT 

r 


/ / 

/ 

sin 

narccos 

0 


\ ' 


sin(0/2) 


\/r'^ cos^{9/2) + oo'^ sin^{9/2) 
It follows from fl7.13p that the kernel of integral eqnation fl7.12p has the form 


(7.13) 


d9. 


Tn{r,u) 


Tn{r,u) 


where the fnnction Tn{r,Lj) is continnous for 0 < a; < r < Ba- Therefore, it follows 
from the theory of Volterra integral eqnations of the second kind that for each a G 
{—B,B) there exists a solntion /3^{oj,a) G C [0,Ba\ of eqnation fl7.12p and this solntion 
is nniqne. Fnrthermore, it is well known from that theory that eqnation fl7.12p can be 
solved iteratively as 


0 d 

^!(^,a) = -— M {^/^{p,a)] {u) 


(7.14) 
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Po 

/ k 1 

'fn \r,a)T„(r,a;)dr - — M {^/J^{p,a)] (u) 


fc = 1,2, 


(7.15) 


and this process converges in the space C [0, pg] to the solntion /9„(a;, a) of eqnation fl7.12p . 
Formnlae fl7.14p and fl7.15p hnish onr second reconstrnction procednre. 
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